Correlated hopping of electrons: Effect on the Brinkman-Rice transition and the 

stability of metallic ferromagnetism 
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We study the Hubbard model with bond-charge interaction ( "correlated hopping" ) in terms of the 
Gutzwiller wave function. We show how to express the Gutzwiller expectation value of the bond- 
charge interaction in terms of the correlated momentum-space occupation. This relation is valid 
in all spatial dimensions. We find that in infinite dimensions, where the Gutzwiller approximation 
becomes exact, the bond-charge interaction lowers the critical Hubbard interaction for the Brinkman- 
Rice metal-insulator transition. The bond-charge interaction also favors ferromagnetic transitions, 
especially if the density of states is not symmetric and has a large spectral weight below the Fermi 
energy. 
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The microscopic origin of magnetic ordering in sys- 
tems like transition metals, transition-metal oxides, and 
high-temperature superconductors is intricate, since it is 
due to correlations between the electrons. The simplest 
model to attempt a description of such systems is the 
single-band Hubbard modeltrEl 



Hubbard 



(1) 



where hats indicate operators. This model describes 
the competition between kinetic and potential energy 
which is at theJieart of the quantum-mechanical corre- 
lation problem.a The Hubbard interaction U represents 
the Coulomb repulsion of electrons in the same orbital 
at a given lattice site. It is given by the matrix element 
U = {ii\V{r — r')\ii) of the Coulomb potential, and is 
typically on the order of a few eV. The matrix elements 
involving neighboring lattice sites i and j are generally 
smaller than U , but may not be negligibly small. One of 
theraisjthe bond-charge interaction Xij = V{r — r') 
\ij)}auu which is typically on the order of 0.1-1 eV and 
hence is comparable in magnitude to the tight-binding 
hopping amplitude i^. It describes a density-dependent 
hopping of the electrons 



(2) 



which only contributes if the lattice site from or onto 
which an electron with spin a is hopping is occupied by 
an electron with spin a ("correlated hopping"). The ef- 
fect of Hx competes strongly with both the kinetic en- 
ergy and the Hubbard interaction. Correlated hopping of 
spin-cr electrons between two sites is enhanced if spin-a 
electrons are present, but this in turn will cost the latter 
kinetic energy, as well as Coulomb energy for the double 
occupations. Moreover, in a band picture the coupling of 
densities and kinetic energy can lead to a band narrowing 
which lowers the amount of energy that is necessary for 
a ferromagnetic spin polarization. These mechanisms ex- 
plain why, in principle, Hx may play an important role in 



the stabilization of ferromagnetism and the localization 
of electrons. 

The model H = i/Hubbaid + Hx cannot be solved ex- 
actly, and only mean=field and finite-size, diagonalization 
results are availableJj Recently, SchillerEI showed how to 
incorporate jEfj^ into the framework of dynamical mean- 
field theory,E3Ea but no numerical results have yet been 
obtained for H . 

One of the standard tools to approach the correlated 
electron problem is the Gutzwiller wave functiond 



1^ 



= n [1 " (1 ~ 5')^«T^a] 1*0 
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where g is a variational parameter (0 < g < 1), and the 
starting wave function |$o) is a product state of spin-up 
and spin-down Fermi seas. By construction both |<I>o) 
and I^'g) are translationally invariant and have a fixed 
particle density n = n| -I- and magnetization m — 
— n|. The uncorrelated case U = Xij — corresponds 
to g = 1 , while U = oo forbids any doubly occupied sites 
and thus corresponds to g = 0. For m ^ and g the 
wave function |5'g) describes an itinerant ferromagnetic 
state. Starting wave functions j^o) with other broken 
symmetries can also be considered. Here, however, we 
will only consider paramagnetism and ferromagnetism. 

Using the Gutzwiller wave function one may, in prin- 
ciple, calculate expectation values of any operator A as 
(A)g = (^'g|^|4'g)/(*g|^'g)- The energy expectation 
value E — {H)q, when optimized with respect to g, is an 
upper bound for the exact ground-state energy of H by 
the variational principle. The variational energy E can 
be written as 

E={H)g^J2 '^"^ka + E 2aRe(a;fe J + Ud, (4) 



where the fc-space occupation n^^, bond-charge occupa- 
tion Xf^^, and double occupation d are defined by (L is 
number of lattice sites) 
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_ J_ \ " ik[R.i-R 
J. 2^ 



(6) 
(7) 



(13) 
(14) 



They each depend on g, n, and m. Here and are 
the Fourier transforms of tij and Xij , respectively. By 
convention, tn = Xn = 0. 

In general the expectation values in Eqs. (||)-(0) are 
not independent of one another. In particular, for the 
Gutzwiller wave function the fc-space and bond-charge 
occupation are closely related. Using the techniques of 
Ref. 13 we obtain, for all Bravais lattices in arbitrary 
dimensions I?, 
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with nl^ 



''kalg- 



^1 as the uncorrelated Fermi function. 
It should be noted that although Xj^^ and ti^^ are linearly 
related, the bond-charge energy and the kinetic energy 
will generally not be, since the regions of the Brillouin 
zone inside (n^^ = 1) and outside (n^^^ = 0) of the Fermi 
sea contribute differently. 

Gutzwiller approximation. — Now we turn to evaluate 
the variational energy E within the Gutzwiller approx- 
imation, which is known to yield the exact evaluation 
of expectation values in terms of the Gutzwiller wave 
functian m. the limit of infinite spatial dimensions {D — > 
oo).tj'll3~tJ It describes a Fermi liquidlj with piecewise 
constant fc-space occupation. 



(9) 



where the discontinuity at the Fermi surface is given by 
[v/("<T - d){l -n + d) + y/{n^ - d) d]^ 



Qa = 



Ua- (1 - ricr) 
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The variational parameter g is related to the double oc- 
cupation d by 



{l-n + d)d ' 



(11) 



and it is convenient to use the latter as variational pa- 
rameter. For the variational energy we obtain 

E = E°- ^(1 - q„){eoa + Co<.) + C/eff (d - n^ni), 

(12) 

where eoo- = (1/-^) X^fe ^fe'^fctr the uncorrelated kinetic 
energy, and similarly ^ocr = (1/^) Z]fc The effec- 

tive Hubbard interaction Ucs and the uncorrelated vari- 
ational energy E° (= i?|g=i) in Eq. ( |l2| ) are given by 



Within the Gutzwiller approximation the bond-charge in- 
teraction thus leaves the form of the variational energy 
unchanged, but enters into the effective kinetic energy 
and effective Hubbard interaction via ^oo- ■ The effect on 
the kinetic energy can be interpreted as a spin-dependent 
band narrowing or widening, which is also found in the 
Hartree-Fock approximation of the correlated hoppiag 
term and can lead to a stabilization of ferromagnetism.u'u 
However the Gutzwiller approximation reveals two dis- 
tinct effects that cannot be resolved in ordinary Hartree- 
Fock theory, where the suppression of double occupancies 
can only be achieved by spin polarizing the system. On 
the one hand, the kinetic energy increases if ^oo- > 0, 
i. e., the effective hopping becomes smaller. This corre- 
sponds to an effective narrowing of the band, which low- 
ers the amount of energy that must be expended for a 
ferromagnetic spin polarization. Furthermore the bond- 
charge interaction contributes to Ucff, which suppresses 
double occupancies already at lower values of the bare 
Hubbard interaction U. The balance between these two 
effects is determined by the optimal variational parame- 
ter d. 

Brinkman-Rice metal-insulator transition. — We now 
discuss the effect of Hx on the Brinkman-Rice transition 
that occutSria the Gutzwiller approximation at half-filling 
{n — l)Jl3'La For convenience we define the strength of 
the bond-charge interaction X by 



(15) 



Note that ^oo- will remain proportional to eoo- according 



to Eq. (|1|) for all densities if tij and Xij have the same 
e. g., if they are nonzero only for nearest-neighbor 



range, 
sites. 



In the limit of D 



both i,i and 



Xij must 



botlLbe scaled as 1/ yZij {Zij is the number of neighbors 
ij),E2l which is compatible with Eq. (p^). The dispersion 
Efc enters only through the density of states (DOS) N{e) 
(which determines eoo-), as expected in dimension D — 
oo. We will consider several densities of states below. 

With the above definition of X, the Gutzwiller approx- 
imation energy for n ~ I and m — simplifies to E' = 
8d(l - 2d){l - X)eo + Ud, where eo = '£a '^Oa < 0. Op- 
timization with respect to d yields a critical value for U, 



[/e(X) = 8|eo|(l-X), 



(16) 



above which the localized state with d = is lowest in 
energy. Hence the Brinkman-Rice transition is moved 
to lower U for X > 0, i. e., the bond-charge interaction 
favors localization. (Only U > and X < 1 will be con- 
sidered from now on.) We find that the U dependence of 
the double occupation d, the discontinuity of the fc-space 
occupation q (= go-), and the energy E is formally the 
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same as in the original Brinkman-Rice theory for X = 0, 
i. e., 
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except that Uc now depends on X [Eq. (|16D]. We note 
that the simultaneous vanishing of q and double occu- 
pation d at a finite value of U is characteristic of the 
Brinkman-Rice transition, in contrast to the numerical 
solution of the Hubb ard m odel, where d remains nonzero 
across the transition.EJIill 

Ferromagnetic transition. — The instability of the par- 
amagnetic state toward ferromagnetism can be deter- 
mined from the bulk susceptibility x- For half-filling, 
we obtain 



1 



X 2iV(ej.) 



{l-X)[l 



pU{U + 2Uc 

{u + UcY 



+ 



rXUc 
U + Uc 



(20) 



where we have introduced the dimensionless parameters 



p = 4iV(ej^)|eo|, 
r = 4iV(eF)ej^. 



(21) 
(22) 



The Fermi energy e^^- in Eq. (£^) represents an absolute 
scale since the first moment of the DOS is fixed at zero 
(due to tii = 0). There are two factors in x that can 
diverge: either g — > 0, i. e., the effective band mass 
m* /m = q^^ diverges at Uc{X), indicating a localization 
transition, or the Stoner-type factor in square brackets 
in Eq. (^0|) vanishes at 



uux) 



y^r^X^ - 4p{l - p){l - Xy 
2(1 -P) 



rX 



X-1, 



(23) 



signaling an instability toward ferromagnetism. The lat- 
ter instability precedes the localization transition when- 
ever p > pfm , where 



Pfn 



1 



rX 



2(1 -X) 



(24) 



These results reduce to the known values pfm = | and 

f/fm = i^p/ip-l) - 1) Uc for X = 000 

Let us first consider the effect of the bond-charge inter- 
action in the case of a symmetric DOS, N{e) = N{—e), 
which results if hopping takes place only between differ- 
ent sublattices of a bipartite lattice. The Fermi energy 
at half-filling is then ep = 0, hence r = 0. In this case 
Pfm = |, the same criterion as for X = 0. On the other 



hand, for an asymmetric DOS ferromagnetism is favored 
by X > if r < 0, i. e., ep < 0. This is the case if the 
Fermi energy is below the center of mass of the DOS, 
which means that there is large spectral weight below 
the Fermi energy. The tendency toward ferromagni 
in such a situation was already proposed long agol 

We now consider nearest-neighbor hopping tij = 
—t*/Vz and bond-charge interaction Xij = X t* j^fZ 
on several infinite-dimensional lattices. Their densi- 
ties of states and phase diagram are shown in Fig. 
The Gaussian DOS of the hypercubic lattice, A^(e) 
= exp(— e^/2)/V27r, has a parameter p = I/tt, and 
DOS of the Bethe lattice, N(t) = 
has p = 32/37r'^, and both have r = due 
Since p < \ \vl both cases. 



the semielliptic 
\/4~^/27r, 



to particle-hole symmetry. Since p < | 
the metal-insulator transition at [4, drawn as a solid 
line in the phase diagram in Fig. |^, will mask the fer- 
romagnetic phase. The variational phase diagram for 
^^Hubbard (i- e., X = 0) oit-thc hypcrcubic lattice was cal- 
culated by Fazekas et al.,Ej who predicted the ferromag- 
netic and antiferromagnetic phases to coexist as the sys- 
tem phase separates and to preempt the metal-insulator 
transition. Here we consider only homogeneous ferro- 
magnetic phases, thus allowing the metal-insulator tran- 
sition to take place, and do not attempt to distinguish 
between paramagnetic and ferromagnetic insulators at 
half-filling, which are degenerate in energy (ii^ = 0). 

A more complicated scenario arises if the lattice sys- 
tem is not particle-hole symmetric, so that the DOS is 
asymmetric and thus r ^ 0. The generalized infinite- 
dimensional fee lattice, with hopping scaled as t,, — 



-1/^2D{D - 1), has a DOS (Ref. 



exp( -(l + V2e)/2) 
7r(l + V2e) 



(25) 



showing a square-root singularity at the lower band-edge. 
The Hubbard model (|l|) on this lattice has been studied 
numeri|Cally by Ulmke within the dynamical mean-field 
theory,Eil who found ferromagnetism at low enough tem- 
peratures and band-filling. For half- filling one has ep = 
-0.3854, p = 1.6157, and r = -1.0272. Thus ferromag- 
netism occurs in this case already for X = 0, but the 
critical U is lowered by the presence of X > 0. Hence a 
ferromagnetic phase is found for Ufm{X) < U < Uc{X). 

Finally, for the class of densities of states with p < 
I ferromagnetism is absent for X = for all U, and is 
only enabled by switching on the bond-chargetiateraction 
X > 0. It is useful to consider a model DOSj23 



Nie) 



1 + VT 



2tt 



(26) 



where the tunable parameter a — —1 ... 1 determines 
the distribution of spectral weight: for a = this DOS 
reduces to the Bethe DOS, whereas for a = 1 there is a 
square-root singularity at the lower band-edge similar to 



3 



the fee DOS. For a = 1 we havep = 1.1353, r = -0.5006, 
while for a = 0.9 the parameters are p = 1.1008 and r 
= —0.2821. Thus in these cases only a metal-insulator 
transition is found for X = 0, but for large enough X a 
ferromagnetic phase is predicted, as shown in Fig. ^ 

Away from half-filling. — Since the metal-insulator 
transition takes place only at half-filling, metallic ferro- 
magnetism occupies a larger part of the phase diagram 
for n ^ I. Figure |^ shows results for the Bethe lattice 
with particle density n = 0.9. Whereas for half-filling the 
strong-coupling phase was insulating, now there is metal- 
lic behavior for all couplings, with ferromagnetism set- 
ting in for large U when X is small, and moving to small 
U when X becomes large. Compared to Hartree-Fock 
theory the Gutzwiller approximation predicts a much re- 
duced region of stability of ferromagnetism, which is due 
to the correlated nature of the Gutzwiller wave function. 
While in Hartree-Fock theory, owing to the lack of corre- 
lations, double occupation can be reduced only through 
a global spin polarization of the system, the Gutzwiller 
wave function describes a paramagnetic state with re- 
duced double occupation controlled by the variational pa- 
rameter g. Furthermore, the Hartree-Fock prediction of 
a phase boundary in a range where either U or X is com- 
parable with the hopping amplitude (see Fig. ||), is not 
consistent with the weak-coupling nature of this approx- 
imation, i. e., self-consistent perturbation theory to first 
order in U and X. Hence we expect that the Gutzwiller 
wave function in general provides a quantitatively better 
estimate than Hartree-Fock theory. 




FIG. 1. Phase diagram for the Hubbard model with 
bond-charge interaction X at half-filling {n — 1). The Brink- 
man- Rice metal- insulator transition takes place for Uc{X) — 
8|eo|(l — (solid line). The dashed lines mark the ferromag- 
netic phase transition for the fee lattice [Eq. (p5|)] and for the 
model DOS of Eq. ( p^ ) for a = 1 and 0.9, respectively. The 
inset shows various densities of states, all with unit variance. 
The lower band-edge has been set to the same value for better 
comparison. 
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FIG. 2. Phase diagram for the Bethe lattice at density 
n = 0.9. The Gutzwiller approximation gives a much smaller 
region of stability of ferromagnetism than Hartree-Fock the- 
ory. 
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FIG. 3. Phase diagram for the model density of states 
[Eq. (p^j for a — 0.9 away from half-filling. 



Figure g| shows the U vs n phase diagram for the model 
DOS of Eq. (|6|) with a = 0.9 for several values of X. Al- 
ready for X — the Gutzwiller theory predicts a large re- 
gion of ferromagnetic ground-atates, in qualitative agree- 
ment with numerical results.E3 The bond-charge interac- 
tion again leads to a further stabilization of ferromag- 
netism. 

Conclusion. — We found that within Gutzwiller's ap- 
proach the bond-charge interaction can enhance the in- 
stability towards ferromagnetism both at and away from 
half-filling. This effect is particularly strong when the 
uncorrelated DOS is asymmetric and there is large spec- 
tral weight below the Ferrijii gnjirgy. This provides further 
support for the conclusionEjIlj that such a situation is fa- 
vorable for ferromagnetism. At half-filling, the presence 
of the bond-charge interaction leads to a metal-insulator 
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transition at lower values than in the standard Brinkman- 
Rice scenario, since it tends to immobilize the electrons. 
Although the Gutzwiller theory can be expected to be 
reliable only at small to intermediate couplings, it rep- 
resents a major improvement over Hartree-Fock theory, 
which for example cannot describe a nonmagnetic local- 
ization transition. 

In conclusion the bond-charge interaction leads to a 
subtle competition between paramagnetism, ferromag- 
netism, and localization. Of course, a variational method 
is not capable of proving the actual stability of a phase. 
It can only provide estimates for the occurrence of insta- 
bilities. Nevertheless, since the Gutzwiller theory treats 
kinetic and interaction effects nonpcrturbatively on the 
same footing, it provides additional insight into the phys- 
ical mechanism behind these instabilities. 
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